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ABSTRACT 


For  an  integer  It  >  1,  let  t  :■  (t. )  be  a  nondecreaaing  real  sequence 
with  t^  <  ,  and  let 


i.*,t 


(x) 


^ti+1'***'ti+k^  “  [ti#...,ti+k_1] )(— x)  + 


k-1 


It  is  well-known  that  N  are  B- splines  of  order  k  for  the  knot 

sequence  t.  Suppose  that  :■  ( )_^  is  a  sequence  of  integers  and 
x  :■  t  .  Then  N.  .  allows  the  following  representation: 

"j.k.t  •  I  ej,k,..tU)»i,k,t  • 

The  coefficient  sequence  P  is  called  a  discrete  B-spline  with  t  and 

J«X»T,t  m 

with  respect  to  t.  This  paper  develops  several  properties  of  discrete 
B-splines  and  proves,  in  particular,  the  total  positivity  of  the  discrete 
spline  collocation  matrix. 


AMS(MOS)  Subject  Classification:  41A15,  15A48. 

Key  Words \  B-splines,  discrete  B-splines,  total  positivity. 
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"Mathematics  Department,  University  of  Wisconsin-Madison,  Madison,  WI  53706. 
Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


SIGNIFICANCE  AND  EXPLANATION 


Spline  interpolation  is  an  important  tool  of  approximation.  Usually,  it 
costs  less  computation  and  yields  a  good  approximation.  The  question  of 
existence  and  uniqueness  of  such  an  interpolant  is  settled  by  the  Schoenberg- 
Whitney  Theorem,  which  is  the  basis  for  spline  interpolation.  There  is  a 
strong  relationship  between  a  spline  and  the  coefficients  in  the  expansion  of 
the  spline  into  a  B-spline  series.  In  many  ways,  the  coefficient  sequence 
behaves  like  the  spline  it  represents.  For  this  reason,  it  is  called  a 
discrete  spline.  In  this  paper,  we  develop  several  properties  of  discrete  B- 
splines  and  prove  the  discrete  analogue  of  the  Schoenberg- Whitney  theorem.  It 
is  expected  that  the  result  obtained  here  would  play  a  role  in  discrete  spline 


interpolation,  discrete  minimisation  and  other  related  areas. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
suonary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


I 


* 
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TOTAL  POSITIVITY  OF  THE  DISCRETE  SPLINE  COLLOCATION  MATRIX 


RONG-QIMG  JIA* 


Discrete  polynomial  splines  on  a  uniform  mesh  were  first  introduced  by 
O.L.  Mangasarian  and  L.L.  Schuoaker  [9] ,  where  discrete  polynomial  splines 
were  defined  as  the  solution  of  certain  discrete  minimization  problems.  Later 
on,  L.L.  Schumaker  [11]  gave  a  description  of  constructive  properties  of  these 
discrete  polynomial  splines.  T.  Lyche  in  his  thesis  [8]  translated  many 
theorems  on  continuous  polynomial  splines  into  their  discrete  analogues. 
However,  they  did  not  view  discrete  B-splines  as  B-spline  coefficients  of 
continuous  splines,  which  allows  consideration  of  discrete  spline  for 
arbitrary  meshes.  It  was  C.  de  Boor  who  first  took  such  a  point  of  view.  In 
my  opinion,  de  Boor's  point  of  view  has  some  advantages  (see  the 
postscript).  Thus  we  shall  develop  de  Boor's  idea  in  this  paper  and,  in 
particular,  prove  the  total  positivity  of  the  discrete  B-spline  collocation 
matrix. 

Let  us  begin  with  some  notations.  As  usual,  X  denotes  the  set  of 
integers,  S  the  set  of  real  numbers,  and  A®  the  set  of  functions  on  B 
into  A.  Thus,  S2  is  the  set  of  real  bi-infinite  sequences.  For  i,j  6  S, 
we  mean  by  [i,j]  the  set  (n  e  Zi  i  <  n  <  j}. 

For  k  6  E,  k  >  1,  let  t  :*  (t^)_—  be  a  non-decreasing  real  sequence 
with  t^  <  t^+lt.  It  is  well-known  that 


Ni,k,t*  "  ( tti+r*#,,ti+k1 


lV,,,#ti+k-11)(*  " 


k-1 
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are  B-splinea  of  ordar  k  for  the  knot  eequence  t.  Here, 


tP0i • • •  *Prl f 


denotes  the  r-th  divided  difference  of  the  function  f  at  the  points 

p  ,...,p  ,  end  (jc-t)*-1  *-  (max  {0,x-t})k 
or  ▼ 

o» 

Suppose  now  that  (u.)_.  is  an  increasing  sequence  of  integers*  For 

j 

T  *■  t  ,  consider  the  B-eplines  associated  with  the  knot  sequence 

j  “J  - 

T  X-  (t.)^  ! 


3  »h,T 


(t)  -  <[T 


j+11 


,Tj+h3  “ 


Since  N.  is  also  a  spline  with  knots  tf  it  can  be  represented  as  a 

j#k,T  ■ 

linear  combination  of  the  N,  's,  by  the  Curry- Schoenberg  theorem  (see 
[4|  113]). 


(1) 


Nj,k,i  ”  [  Sj#K,T,t(i)  Hi,k,t  * 


Following  de  Boor  [2],  we  make  the  following  definitions 

Definition  1.  The  coefficient  sequence  B,  ,  _  x  €  R*  in  (1)  is  called  a 
discrete  B-spllne  with  knots  t  and  with  respect  to  t. 


It  is  known  from  [2]  that 


<»  •j.k.i.t'11 


]  (  •  •  t  ) 
j+kJ  1 


*  ti+k-1)+ 


r 


When  k«1,  (2)  reads 


-2- 


# 


(2*) 


"j.i.j.s'1' 


(t 


j+1 


V+ 


where 


(t4- 


V+ 


1  if  t  >  i, 
0  if  t  <  i. 


If  t  and  t  are  clear  from  the  context,  0  will  be  abbreviated  as 

*  *  J 

0  ,  or  even  to  0.  . 

j  *R  3 

Remark  1.  Definition  1  uses  a  different  normalization  than  do  (5.10a) 


and  (5.10b)  of  [2].  Clearly, 


W*U) 


i+k 


- 1,  V1,# 


where  a^d)  is  in  the  sense  of  (5.10b)  of  [2]. 

Let  us  now  establish  some  basic  properties  of  discrete  B~ splines. 
Lemma  1.  (Narsden's  Identity) 


(3) 


Proof.  It  follows  from  Marsden's  Identity  (see  [21)  for  continuous 
B-splines  and  (1)  that 


\  "i,k,$  “  1  “  |  *3,k,T  “  Bj,k,T,|(i)Hi,)Crt^ 


1  -  I  N 
i  3 


Since  R  (i  e  S)  are  linearly  independent,  (3)  must  hold. 

1#K,T 

Lemma  2.  (Composition  formula)  Let  t  »•  (t^)  n  be  a  non-decreasing 
real  sequence  with  t^  <  t^+Jt<  g  a  subsequence  of  t  and  t  a  subsequence 
of  £•  Then 


-3- 


(4) 


.(8)8,  u  A  k  • 

j*k»T,t  *#*#£#5 


Proof .  By  ( 1 ) , 


I  “j.k.i.t111  "l.k.t  ■  "i.k.j  ■  £  Bj,k.:,eu>"t.k,e  ■ 


Since  I.  ,  .  (i  «  I)  are  linearly  independent,  (4)  follows  fro*  the 
*#*#5 

above  equality. 

3.  (Recurrence  relation)  For  k  >  2  , 


*1+1  k-1(i)  k-1(1) 

(5)  B.  „(i)  -  <T__  -  t__  J  - +  <t__  -  Ta>  r-3tiL- 


i.*' 


j+k  -i+k-1'  tj+k  -  Tj+1  '  '“i+k-1  i'  xj+k_,  -  tj 


Proof.  Mote  that  for  any  T  €  t. 


11  '  Vi>.  — IT  -  Vk-2>»  (T  -  Wi>* 


t(T  *  ti-f1,+  ***<T  "  ti+k-2>+J(T  ~  ti+k-1> 


Applying  Leibniz’ •  formula  to  the  above  product,  we  obtain 

*j,k(i)  *  (Tj+k"Tj)lTj#***#Tj+kI(*"ti+1,+  ••*<*~ti'Hc-1,'»> 

J?k 

*  . V'"*!.!1.  ,"tl.k-2>.,I,r . 

(Tj+k"Tj,lTJ,'"'1j+ltI  ,t'’tl+l,+  " ' ( ”tl*k-2,+<Tj>k"ti^k-1 1 

»  . Vk-l""W* 


(  X  j+k“*i+k-1  >  1 1 T  j+1  *  •  •  •  '  T  j+k5  “  1 X  j  *  •  •  •  •  T j+k-1 1  * (  *+  *  *  *  (  *"*1+^2  >+ 

+  <‘lj+k"Tj,Itj'***»Tj+k-1I(*"ti+1>+  *"ti+k-2,+ 

t(T  j+k"Vk-i  > c  Vi . V*1  +  ( WrTj  > 1 V — '  Tj+k-i J }( -Vi  K— <  -  Vk 

Vi  k-i(i)  \  k-i(i) 

(t  -  t  )-3ridri —  +  /t  -T  )  — ^K  - —  . 

,Ti+k  li+H  t  -  T  tx:i+k-1  v  T  -t 

3  *  1TJS  1  l^+K  ‘  j+1  j  1  j+k-1  ‘j 


This  proves  Leona  3. 

4.  For  a  fixed  s  6  I,  let 


v»  *" 


t  if  t  <  ■ 
*  +  1  if  t  >  8. 


>*  tv  a&d  g  j«  •  In  other  words,  g  is  formed  by  dropping  an 

entry  fron  t.  Then 


0,  .  .(i)  ■  0  for  £  <  i-1  or  t  >  i> 

t.k.g,t 


(6a) 

(6b)  6.  .  .  . (i)  >  0  with  strict  ineauality  iff  t._.  >  t  » 

i-1,K,g,tf  11  "  .  ■  i  -  .  i  ■  1 1  ■■  i+k  e 

(6c)  0.  k  o  t*1^  >  0  with  strict  inequality  iff  t .  <  t  . 

9  9m9m 

Proof#  If  tt  >  t  0  then  H§  .  ^  ,  .  /  and  it  follows  that 

a  s  %#K|P  *▼ 1 1 k  f  t 


et#k,g#S(i)  “  6i,Jt+1  for  a11  i#** 


in  the  sane  way,  for  *Vk+1  <  t#  , 


*»  v  «  *.(i)  “  *11  i,*. 


Mow  (6a)  is  easily  derived  fron  what  has  been  proved.  Moreover,  when 
*t  <  *.  ‘  Wi-  *• h*” 


S' 


Hi,k,g  *  *i,k,£,t(t)Ht,lc,t  +  ®t,k,g,tU+1)  Hi+1,k,t  * 

It  is  known  that  fc  has  the  same  sign  as  ffjj  ^  ^fcU)  in  (tt#  t^  +  e) 
for  sufficiently  snail  e  >  0  (see [2]),  so 


similarly 


®l,k,g,t(t)  >  °  £°r  <  *.  <  Vk+1, 


*t,k,e,t  (1+1)  >  0  £or  **<*.<  Wi  • 


nil— hi  1  rlnj  these  facts  we  get  (6b)  and  (6c). 


5.  Suppose  £  €  S*  and  x^  -  t^  for  all  j  €  S.  Then 


‘j.x.I.S11’  '  0 


with  equality  if  and  only  if  one  of  the  following  four  cases  occurs » 


(7a) 

ti< 

t  > 

**3 

j 

(7b) 

ti“ 

t  and 

Wj 

(7c) 

*i+k 

>  \ 

Mj+k 

(7d) 

fci+k 

“  % 

H  j+k 

t.  -  t  and  max  {p|  t  -t. }  >  max  {q|  t  «t  }| 
i  Wj  i+P  i  Pj+q  Pj 


and  nax{p|  t-t}  >  nax{q|  t  «t  } 

p  j+h-q  j+k 


Proof*  We  use  the  linear  functional  X^  given  by  the  rule* 


X  f  «-  l  (-1) 
1  r-0 


k-1-r  T(k-1-r)(5)Drf(?)^  aU  f# 


gSjgwr  <*v 

■ 


By  the 


where  f(t)  ***  (ti+1  “  t)...(ti+J{-1  ~  t)/(fc-1)l  end  <  £  <  ti+k  . 
de  Boor-Fix  Theorem  (see [4;  116-118}), 


(9) 


Xi*j,k,j  *  6j,k,T,tU)’ 


In  case  (7a),  choose  £  so  that  t.  <  £  <  t  .  Then 


N.  k  t(5)  "  K  v  r(V  “  •••  *  M5lt*1i($)  *  °' 


Hence  (8)  and  (9)  yield 


“  XiHj,k,t  "  °* 


In  case  (7b)  write  c  *«  max  {p|t..  •  t. ),  d  »■  max  {q|t  “  t.  }.  Then 

"*  1  “)♦»  "l 

d  <  c— 1  and 


*Wv>  -  "i.k.j'v  -  -  ■  "m£"«v>  - 0 

-  ...  -  t^0”1  * (tj+)  -  o. 


Taking  £  -  t^+  in  (8)  and  substituting  these  values  into  (8),  we  obtain 


'j.k.T.j'1’  '  Vj.k.J  ■  "• 


Cases  (7c)  and  (7d)  can  be  treated  in  the  same  way. 

Now  suppose  that  none  of  (7a)  -  (7d)  is  true.  We  want  to  show 
»4  v  _  Ji)  >  0.  Let 


-7- 


B  im  {t|  Vi j  <  l  <  le  jjl  and  |B|  i*  the  cardinality  of  E. 


We  shall  proceed  by  induction  of  |e|.  The  case  |s|  ■»  0  is  trivial.  The 

case  |E|  «  1  is  reduced  to  Leona  4.  Assume  now  that  our  statement  is  true 

for  |E|  <  n.  We  want  to  prove  our  statement  is  also  true  for  |E|  “  n.  Take 

any  s  e  E.  Let  g  be  defined  as  in  Lemma  4*  that  is,  -  t  for  l  <  s, 

v  ■  i+1  for  1  >  s  and  p.tm  t  .  By  Lemma  2  and  Lemma  4, 

1  1 


B.  .  r  Ai.)  m  l  6,  v  T  ft<*>B,  .(i) 


B,  w  k  a  „(i>  +  0,  k  T  (i) B.  v  A  .  (i). 

j*h,T,g  i  1,k,g,t  j#k,T,g  i#k,g,t 


All  terms  that  appear  in  the  above  equality  are  nonnegative.  It  seems 

appropriate  to  treat  the  following  three  possible  subcases  individually. 

(i)  t^  >  t^  .  In  this  case,  »k,g,t(i)  >  0  by  (6b).  We  need  to 

show  B.  .  ,  (i-1)  >  0.  If  i-1  >  s,  then  v.  .  -  i  and  v.  ...  «  i  +  k; 
3,k,t,g  i-i  i-i+k 

so  B.  .  _  (i-1 )  >  0  by  induction  hypothesis.  Assume  now  i-1  <  s.  If 

t  >  t  ,  then  t  >  t  or 

j  j 

fci-1  “  t|iJ  *nd  “*  ^'Vi  “  Vi1  "  0  <  “*  {qltMj+q  ’  tMJ}  * 

Hence  B.  .  (i-1 )  >  0  by  induction  hypothesis  again.  Finally,  suppose 

J,K,i,g 

t.  *  t  .  Then  t.  >  t  >  t  implies  t  -  t.  . 
i  “j  i  «  *  «  i 

otherwise  >  i  and  i  <  s  would  imply 


max  (plt1+p 


t  >  >  max  (q| t 

j+q 


a  contradiction.  In  conclusion, 


Thus  <  if  for 


■**  lpl»i.,+p  -  »i-,>  -  *“  w*n,  '  'i*  ‘  *“  fa'Sj  ■ 


so  that  8.  .  n(i-1)  >  °* 

j|K|T|g 

(ii)  t  <  t  .  This  case  can  be  treated  in  the  same  vray  as  (i)  is. 

'  i+k  s 

(iii)  tA  <  tg  <  ti+Jc  .  Lemma  4  tells  us  that  both  ,k,g, tll) 

0  (i)  are  positive  in  this  case.  Thus  we  need  to  show  that  at  least 

i»k,£,£ 

one  of  B ,  1 4.-1  >  and  84  .  ,  (i>  is  positive.  If  either  t.  >  t  or 

j#k,t,g  3#h,T,g  1  Vj 

t  <  t  ,  then  this  holds  by  the  observation  made  in  (i).  Next,  suppose 
i+k  Mj+k 

*»•  ■  ‘B)  ■  v*  -  v  - 


either  max  {p|ti+p  -  t^  <  max  {q|ty  *  t^} 


or  max  (pit.  .  «  t  }  <  max  {q|t  -  t  }, 

1+K  p  1  yj+k-q  j+k 


then  one  can  easily  get  0.  .  (i-D  >  0  or  B ,  .  D)  >  °»  using  the 

JrKfTfg 

same  argument  as  in  (i)e  The  remaining  case  to  be  discussed  is 


{plti  -  tL)  *  max  {qlt^  =  t±>  and 


max  {plti+k.p  -  ti+k)  -  max  tqlt^+^  -  t^}. 


Let  c  :»  max  {p|ti+£>  -  t^,  e  :  =  max  tplti+k_q  =  ti+lt^  Then 


hence 


Vc+1  >  i+c.  Uj+k_e_,  <  i+k-e. 


j.  .  “  <  ( i+k-e- 1 )  -  (i+c+1)  *  ( j+k-e-1 )  -  (j+c+1). 

j+k-e-i  3+c+i 


This  means  «  e  jj(  which  contradicts  the  choice  of  s.  Lemma  5  is  proved. 


4 


We  are  now  in  a  position  to  prove  our  main  result. 

Theorem  1.  Let  t  *■  be  a  nondecreasing  real  sequence  with 

t^  <  all  It,  an  increasing  integer  sequence,  »-  ty  and 

MM  ^ 

let  x  »■  (t^)^  •  Let  be  the  sequence  of  discrete  B-splines  of 

order  It  with  the  knot  sequence  T  and  with  respect  to  Let 

i.  <  i„  <...<  i 
12  m 

be  a  finite  increasing  subsequence  of  integers,  and  set 

0,-  <urj>  -CBjCV),  <r<B. 

Then  for  every  subsequence  q^  <...<  q^, 

-  s 

LV  •••» 


(10) 


det  0 


i^  ,  * « *  i 


>  0 


with  strict  Inequality  iff  both  of  the  following  conditions  are  satisfied; 

(i) .  0  (i  )  >  0  for  all  r  -  1,2,...,m? 

qr  r 

(ii) .  If  there  is  soae  s  e  such  that  tA  *  t^  for  some  r,  then 


where 


Proof .  Write 


Ar-d  <  Ar  -  dr 
r 


lr  t-  k  -  max  {p|ti  -  t±  > 


+p 

r  r 


A  »«  0 


i1 ' *  *  * '  Sn 

ql'*,,,qm 


If  0  (i  )  -  0  for  some  r,  then  0  (i.)  *  0  for  all  A,  j,  with 

qr  qj  * 

1<£<r<j<m,  by  Lemma  5.  Thus  columns  r,...,m  of  A  are  linearly 


dependent  and  det  A  *  0.  Without  loss  of  generality  we  may  assume  further 


that  both  the  first  superdiagonal  and  subdiagonal  of  A  are  positive. 


10- 


(11)  0  (i...,)  >0,  r  ■  1,...,m-1  and  6  (1 _ .)  >  0,  r-2,...,m. 

q^.  r+i  q^.  r-i 


Otherwise,  we  would  have,  aay,  0  (i  )  ■  0  for 

qr+1 ,1c  r 

that  0  (i^)  -  0  for  any  i,  j  with  1  <t<r<j<m.  Thus 


r.  It  would  follow 


i^  . * • » »ij 

ir+r***,im 

det  A  -  det  U 

•  det  0 

q1  #•»» 

L^i . \ 

where  det 


i1#*««*ir“l  t  tf  rir+1,*,,,i»l 
Lv-'vJ  L’v+i . sJ 


are  lower  order 


L*r+1 . V 

determinanta  of  the  same  fora.  If  a  »  1,  then  X  is  a  1x1  matrix  and 


det  A  >  0,  trivially.  Thus  if  we  use  induction  on  m,  then  det  A  would 
already  have  the  property  declared  in  Theorem  1.  Proa  now  on  we  always  assume 
(11)  to  hold. 

We  point  out  that  (11)  yields 


(12) 


0<J^(ir±1>  >0,  r  -  1,...,m. 


The  only  thing  we  have  to  prove  is  0  (i  -1)  >  0  and  0  (i  +1)  >  0,  while 

*1  1  ^m  * 

in  all  the  other  cases  this  is  a  direct  consequence  of  (11).  Since 


0  (i.)  >  0,  i.  >  u  ,  hence  i  -  1  >  u  -  1  >  |i  ,  then 

q2  1  1  q2  1  q2  q1 
u  lr(ii  -  1)  >  o,  obviously.  If  t  *  -  t  and  t,  ,  <  t.  , 

V*  1  l1  1  V  41  1  A1 


then 


max  {pl^  _1  *  t£  -1+p}  -  0  <  max  {q|t 


.  *  V1’’ 


So  we  also  have  0  (i  -1 )  >  0.  The  last  possible  case  is 
q1  1 

t.  ,  *  t  ■  t.  .  Then 

V’  ‘i 


«ax  {p|t 

1 


—  •  v-Ul,>  ■ 1  *  ■**  (plti, . *iV 


<  1  +  max  {q|ty  -  ...  -  ty  >  <  max  {qjt^  •  ...  -  t^  )i 


q2+q 


q.+q 


-ii 


r 


therefor*  0  (1,-1)  >  0.  Similarly,  0  (i  +1)  >  0. 

q1  1  * 


As  In 


5,  let 


St  -  Ul  uq  <  t  <  W  fc,  1  *  H* 


We  will  proceed  by  induction  on  |8|.  If  |E|  “  0,  then  A  is  a  diagonal 
matrix,  so  the  proof  is  trivial.  Suppose  now  that,  for  |B|  <  n,  our  theorem 
is  proved,  and  we  want  to  show  the  conclusion  of  Theorem  1  also  holds  for 
1*1  -  »• 

We  have  proved  that  if  (i)  is  violated,  then  det  A  *  0,  Suppose  now 


that  (ii)  does  not  hold.  Then  there  is 


set  such  that  t.  ■  t  and 

-  is 

r 

4.  Let 


i _ .  “  i  -d  .  Fora  o  by  dropping  s  from  t  as  w*  did  in 

r-a  x  r  •  * 

r 


V  **  ^8i,k,o,t(ir^  1<r<m 


u  4*<|i 


Vk-1 


1<r<m 


Them  A  -  VW  by  Imm  2.  Further,  the  Cauchy-Binet  formula  (see[6l)  gives 


(13)  det  A  »  l  det  V 


5i<?2<..*<{. 


W 

5l,52' 


det  W 


W 


Since  t  -  t  and  t.  “  t.  -d  ,  we  have 

r  8  r-dg+k  r 


8i  k  o  k11,*  "  0  *nd  -a  k  «  >  “  0 

r*k»fi»£  r  Ar "*r “1»k#£*S  *“dr 


t»y 


4.  Furthermore, 


* 


12 


(V  -  0  £or  i  <  W1  or  i  >  4r 


i 


h  •  i  — d  *  •  •  •  i  i  « 
r  r  r 


Consider  the  following  matrix  with  d^+1  rows: 

■'  8ir-dr-1#)c,g,t^ir-dr^  8ir-dr,k/g,t^ir-dr^***6ir-1,k,g,t^ir-dr^  ^  ,k,£,  t  ^r-d^.  ^ 


•  B  (i  )  6  (i  )  •••£  (1  )  B  (t  ) 

Pir-dr-t,k,g,t'  r'  pir-dr,k,fi,tl  r;  pir-1 ,k.£,t' V  pir,k,£,tlV 


All  its  entries  except  those  in  columns  i^-d^ , , . . ,1^-1  ere  zero.  Thus  the 
rank  of  this  matrix  is  no  bigger  than  d£.  Hence  the  d£+1  rows  of  this 
matrix  are  linearly  dependent.  This  shows  that  the  rows  r-d^. ,  r-dr+1,...,r 
of  the  matrix  A  are  linearly  dependent.  Thus  the  rows  r-d£,  r»dr+1#...,r 


of  each  V 


rS'S . lml  . 

M2 . «• 1 

A1  fi2 . Sa] 

“ v  L*.* . <J 


are  linearly  dependent,  so  that 


0  for  all  £,  <  C,  <  ...  <  5  • 

12  ni 


Therefore  det  A  *  0  by  (13). 

Suppose  now  that  both  the  conditions  (i)  and  (ii)  are  satisfied.  We  want 

to  show  det  A  >  0.  We  shall  argue  by  induction  on  |E|  again.  Take 

s  e  E*  Form  v  and  g  as  we  did  in  Lemma  4.  Let  V  and  w  have  the  same 

meaning  as  above.  By  induction  hypothesis  and  Lemma  4,  all  products  that 

appear  on  the  right-hand  side  of  (13)  are  nonnegative.  Let  r  be  the  least 

integer  such  that  t.  >  t  .  Then  t.  <  t  .  There  are  two  possibilities 

Li  8  Vi  * 

to  be  discussed: 


i  ,  <  i  -1 
r-1  r 


In  this  cans,  we  choose 


Si  *" 


i^-1  for  h  >  r, 
i.  for  h  <  r. 


Than  C,  <  5,  <  ...  <  {  .  By  Lawn  4  and  tha  choica  of  tha  ('i, 

1  «  B 


ith.v,t.SUi>’  *  h  •  1,2 . 


In  addition,  if  h  <  r,  then  we  have  v.  -  i.  and  v_  •  i  or 

Si  h  Si+k  *** 


ih+h+1*  Thu*  M2)  together  with  Laoma  5  tella  us  that 


8  k  _  >  0. 

Similarly,  if  h  >  r,  then  we  have  •  i^  -  1  and  »  i^)fc  -  1  or 

ih+k  so  the  above  inequality  also  holds.  By  induction  hypothesis  we  assert 


det  V 


VV'*"Sa 

.  ,  -  >  0  and  det  W  >0. 

*1,*2***,#Ss  r<*2' *  ’  * '**« 

rvs . i.i 

re  have  det  U  >  0. 

ISj  ’  *  *  'Si  J 


By  (13)  and  (14)  we  have  det 


(0)i  -I  *  i  4#  In  this  case,  condition  (ii)  gives 
r  r-i 


i  .  <  i  -  d 

r-d  r  r 
r 


where  d_  i*  )c  -  max  {pit. .  ■  t,  }•  There  exists  an  integer  c, 

r  i+p 

1  <  c  <  d  ,  such  that  i  ,  «.  ■  i  “  (c-1)  but  i  •<  i  -  c.  Thus 

r  r-tc-i j  r  r~c  r 


Vc  <  Ar  -  c  *  (Vc+,}  '  1  '  Wl  -  ’* 


1 


Let 


(16) 


for  h  >  r-c+1 , 
for  h  <  r-c. 


Proa  (15)  and  (16)  we  see  that  5^  <  5  <  <  €a«  Ho*  Le— a  4  yields  that 


Using  the  saae  argument  as  in  (a),  we  get 


det  W 


>  0. 


This  proves  our  theorem. 

Remark  2.  if  t  is  a  strictly  increasing  sequence, 
stated  as  follows: 


•ja.j.s'1' 


>  0 


» 


5  can  be 


with  strict  inequality  iff  t.  <  t . ,  t...  <  t . Furthermore,  in  Theorem  1, 

J  1  X+K  J+K e 

the  condition  (ii)  is  automatically  fullfilled  as  long  as  (i)  holds. 

Remark  3.  it  is  interesting  that  the  Schoenberg-  Whitney  Theorem  (see [9]) 

can  be  derived  from  our  Theorem  1.  Indeed,  let  t  *  (T. )  "  be  a  non* 

■  3  "" 

decreasing  knot  sequence ,  <  • • •  <  T  ,  and  let  t  be  a  knot  sequence 

formed  by  adding  some  knots  to  t  so  that  t  has  exactly  k  multiples  at 
each  T^,  i»1,2,...,m.  Then 


Nj#k,J<Ti)  "  aj,k,T,J(Ti) 


according  to  (1).  Wow  one  could  easily  see  that  the  Shoenberg- Whitney  Theorem 
is  a  consequence  of  Theorem  1. 
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Postscript 

This  work  was  dona  in  July  1980.  Later  X  became  aware  of  the  three 
related  papers  11],  IS]  and  [7] .  Essentially  whether  explicitly  or 
implicitly,  these  three  papers  view  discrete  B-*splines  as  the  coefficient 
sequences  associated  with  the  expansion  of  continuous  polynomial  splines  in 
B-splines.  This  is  just  de  Boor's  point  of  view  (see [21 ).  In  [1],  the  author 
provided  an  algorithm  for  further  subdivision  of  a  knot  sequence.  The  basic 
idea  of  [1]  is  to  investigate  what  happens  when  one  inserts  new  knots  into  a 
given  knot  sequence.  The  essential  idea  of  the  present  paper  is  also 
"inserting  new  knots”  and  "inserting  one  new  knot  each  time".  In  [5],  the 


authors  develop  more  properties  of  discrete  splines. 


1  and  3,  a  part  of 


Lemma  4  and  5  of  this  paper  overlap  with  [5] .  However,  [5]  is  based  on  the 
recurrence  formula,  while  ay  Theorem  1  does  not  need  recurrence  formula  though 


the  proof  for  the  recurrence  formula  (Lemma  3)  is  more  straightforward  in  my 
opinion.  In  [7],  the  authors  give  the  shortest  way  to  prove  the  variation 
diminishing  property  of  B-spline  approximation  by  using  a  geometric 
observation.  Their  methods  can  be  easily  carried  to  proving  that  the 
associated  discrete  spline  collocation  matrix  is  sign  regular,  but  it  seems 
hard  to  detemine  which  minor  is  really  positive  along  this  way.  In  the 


present  paper,  by  the  composition  formular  ( 


2)  and  the  Cauchy-Binet 


formula,  we  are  able  to  obtain  the  exact  criterion  for  the  positivity  of  a 
given  minor.  I  believe  that  the  determination  of  such  positivity  is 
significant  and  expect  that  Theorem  1  will  play  a  role  in  discrete  spline 
interpolation,  discrete  minimisation  and  other  related  topics. 
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B-splines  and  proves,  in  particular,  the  total  positivity  of  the  discrete 
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